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We consider a cylindrical metallic magnet that is set into rotation about a horizontal axis 
by a falling mass. In such a system the magnetic field will cause a radial current which is 
non-solenoidal. This leads to charge accumulation and a partial attenuation of radial current. 
The magnetic field acting on the radially flowing current slows down the acceleration. Using 
newtonian dynamics we evaluate the angular velocity and displacement. We explicitly show 
that the electromagnetic angular momentum must be taken into consideration in order to 
account for the change in angular momentum due to the external torque on the system. 
Further the loss in potential energy of the falling mass can be accounted for only after 
taking into consideration the electrostatic energy and the Joule loss. We suggest that this 
example will be of pedagogical value to intermediate physics students. A version of this is 
scheduled to appear in American Journal of Physics 2011. 
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INTRODUCTION 



Introductory and intermediate treatments of electromagnetism usually shy away from a formal 



m 



Discussions if any are kept at 



discussion of electromagnetic angular momentum and energy. [ 
a qualitative level. Specific examples involving forces and torques as well as the associated work 
done and the change in angular momentum are conspicuous by their absence. In a series of articles 
Gauthier has attempted to redress this lacuna. P 



m 



In one of his articles Gauthier considers a metallic cylindrical shell free to rotate about the 
horizontal axis. [J] A string wound around the shell with a mass attached to its free end causes it 
to rotate with acceleration. A charge is smeared on the shell making it analogous to a solenoid 
carrying an increasing time dependent current and associated magnetic flux. One can show from 
Faraday's law that a retarding torque will act on the shell. This torque is of purely electromagnetic 
origin. Gauthier goes on to demonstrate in a pedagogically satisfying fashion that the system's 
magnetic energy and electromagnetic momentum are stored at the expense of their mechanical 
counterparts. 

In the present article we consider a cylindrical magnet free to rotate about the horizontal axis. 
We treat the magnetic field of the magnet to be constant and pointing along the axis. A string 
wound around the shell with a mass attached to its free end causes it to rotate with acceleration. 
We assume that the magnet is a good conductor and hence a radial current will be set up. The 
magnetic field will act on the radial current giving rise to a counter torque which will oppose the 
torque due to gravity. [6] The problem is complementary to Gauthier's. Instead of Faraday's law we 
need to invoke the Lorentz force, Gauss's law and the equation of continuity. But the perspective 
is the same as Gauthier's: an understanding of the problem involves a detailed consideration of the 
electromagnetic energy and angular momentum. We thus present another pedagogical example for 
an intermediate electromagnetic course. 



II. DYNAMICS OF A ROTATING MAGNET 

Consider a cylinder with moment of inertia I, radius R and length /. Let mass m be hung by 
a light thread wound on the surface of the cylinder. We adopt the cylindrical coordinate system 
{r, (p, z} and unit vectors {f, (j), k} with k being the axis around which cylinder is free to rotate. 
Figure 1(a) illustrates the setup and Fig. 1(b) the cylindrical coordinate system we have adopted. 
When mass m is released, the cylinder rotates with instantaneous angular velocity uj'{t)k. We 
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(a) (b) 
FIG. 1. (a) Setup, (b) The cylindrical coordinate system. See text for details. 



assume no slip condition v = roj' . In that case Newton's second law for rotational dynamics can 
be used to show that angular velocity of the cylinder is 



(1) 



where we have introduced a dimensionless constant 



I + m,i?2 ■ 



(2) 



and a time constant 



9_ 
R' 



(3) 



We now replace the cylinder by a metallic cylindrical magnet in the setup. Let the conductivity 
of the cylinder be a and the magnetic field be uniform and Bk.[6] This is similar to the setup 
described by Gauthier [J] except that we have a solid cylindrical magnet instead of a non-magnetic 
charged cylindrical shell. As the magnet rotates due to the gravitational torque of the falling mass, 
the free charge in the rotating magnet experiences a radially outward Lorentz force. The magnetic 
field will act on this radial current and will produce a counter torque (rem) which will oppose the 
torque due to gravity. The subscript indicates that the torque is of purely electromagnetic origin. 
We can write the equation of motion for the system as 



mgR — Tem{t) = {I + "ii? 



2\ 



dt 



(4) 



In what follows we will determine rem and solve the equation of motion. 

Consider a free charge q at a distance r {r < R) from the axis of the cylinder rotating with 
angular velocity uj. It is subjected to a radial Lorentz iorce our Br. Recall that in literature 
the effect of this force is discussed in terms of motional emf. [7|, [S] From Ohm's law this motion 
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would lead to a current density Jm = aojrBr. It is important to note that this current density is 
non-solenoidal 

r ar 

and hence charge must accumulate. This accumulating charge will give rise to a electric field (Ea). 
Let p{f, t) be the charge distribution due to accumulation. Associated with this p{f, t) is an electric 
ield which can be determined by Gauss's law with the standard pill box construction in the magnet 



1 r 

Ea = / p{r',t)r'dr' (r < R), 

^or Jo 



where we assume that the length of magnet / ^ R. Once again we use Ohm's law to obtain 
Ja = crEa'r. The resulting current density across the magnet is then 

J = a(ujrB^ / p{r' ,t)r'dr'j r, (5) 



^or Jo 

It consists of two terms: the first is due to the motional emf and second is due to charge accumu- 
lation. Next examine the divergence of current density 

V.J = -^(r J) = a (lioB + —p{r, tU 
r ar \ eo / 

From equation of continuity 
we obtain 

I + ^ = -2Ba.it). (6) 
ot eo 

We assume that initially the cylinder is neutral {p{r,t = 0) = 0). Using this, Eq. ([6]) can be 
solved to yield 

p{r,t) = -2Bae-^/^^ [ u{t')e^' /^^dt\ (7) 

JO 

where T2 = eo/cr. The reader may have encountered T2 earlier. When a lump of charge is deposited 
on a conductor, it spreads out in a very short time and with time constant T2. Note that p 
is independent of r. In order to maintain electrical neutrality there will be an accumulation of 
charges at the surface with the sign opposite to p. Hence 

Ea = Q {r>R). (8) 
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Consider an annular element of radius r < R m the cylinder (see Fig. 1 (b) ) . The force {dF) on 
the annular element can be written as 

dF = 2TTlrJdr x B 
= 2TTlrJB{-4>)dr 

The associated torque (dTem) on the ring will be 

dTem{t) = fxdF = 2TTr'^lJB{-k)dr 

The subscript indicates that torque is of electromagnetic origin. The net torque is obtained by 
straight forward integration 

dTem{t) = fxdF = 27rr^lJB{-k)dr. 

Using the expressions for J (Eq. ([5])) and p (Eq. ([7])) one may obtain the net torque after straight 
forward integration 

remit) = ^-^^^ {u:{t) - le-*/- j\{tY'/--dt^ . (9) 

We pause to note that since the charge density is independent of r, the electric field from Gauss's 
law will have a linear dependence on r. Ohm's law implies that the current density would also be 
linearly dependent on r. It now remains to determine the time dependence. Differentiating above 
Eq. (jl]) with respect to t gives 

d^uj 1 d 



Using Eq. © 

d irlaB'^R^ / du 1/1 



(10) 



dt^""" ~ 2 



f--e-*/- / u;{ty/^-dt' + u;{t))) (11) 

dt T2 \ T2 Jo J J 



The integral term on the r.h.s of Eq. (jlip can be eliminated using Eqs. ([9]) and Hence we 
may rewrite Eq. (jlOp as 

d'^uj f TilaB'^R^ l\duj_ mgR 

'd^^\2{I + niR'^) ~dt ~ T2{I + mR^) ' ^ ^ 

We define a time constant 

1 7r/BV/?4 



rs 2(1 + mR'^) ' 



(13) 
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and an effective time constant 

1 _ 1 1 

Te T2 Tg ' 

One can solve Eq. (|12p using initial condition u:{t = 0) = and duj/dt = j^/rf to obtain 

uit) = 4* - (14) 

where f{t) = t-Te{l- 6"*/""=) (15) 
Note that the first term in Eq. (jl4p is same as Eq. ([2]). Also 

f{t) > 

This indicates that the angular velocity in case of magnet is never greater than in the case of 
non- magnetic solid cylinder. We consider two limits of the Eq. ([1] 



Uj{t » Te) W (17) 

< 4* = ^' (18) 

This shows that the effect of the metallic magnet is to decrease the angular velocity from the 
non- magnetic value {lo'). A straightforward integration also yields the angular distance traversed 



T{T3 2rfr2 



^(t) = zrtfii) + ijt^ (19) 

Once again we note the limiting cases 



</'(* < ^e) = 



rfra 2rfr2 
^ 2rr 

which indicates that the magnetic field slows down the rotation. Using Eq. (|14p in Eq. ([7]) and 
integrating we obtain charge density 

pit) = =^^m. (20) 



We remind the reader that f{t) > and hence charge density is negative. Also, as stated earher, 
a compensating positive charge cohects on the surface so that the magnet is on the whole neutral. 
We have already pointed out that p is independent of r and therefore the electric field is 

^. = f^f = -^/Wr {r<R). (21) 
= (r > R). (22) 



We then calculate the total current density 



J=^{t-f{t))r. (23) 

while the net electromagnetic torque on the system is 

rem = ^!-^ -{t-f{t)). (24) 

We pause to comment briefly on the time dependence of above quantities. Note that t — f{t) = Te(l - 
As physically expected both the current density and electromagnetic torque are initially zero and 
they both saturate to a constant value at large times. 

III. VERIFICATION OF ANGULAR MOMENTUM 



The electromagnetic angular momentum (Lgm) is defined in terms of electric and magnetic 



field. 



7|] In our case we have 



Lem = / eor' X {Ea X B) Inr'ldr' (25) 
Jo 



Using the expression for electric field due to the charge accumulation (Eq. (|2T]) ) we have, 
l^em = / t / r'f X r'f X kydr' 

T\^T2 Jo 

T-r-r2 Jo 

= f[t)k (26) 

Since f{t) is always positive it means that some angular momentum is transferred to the elec- 
tromagnetic field and consequently the mechanical angular momentum is diminished. One can 
check that differentiating the above equation for Lem, we will obtain precisely Eq. (124p for elec- 
tromagnetic torque. The mechanical angular momentum of the system is 

L„ = (I + mR^)ujk = {I + mR^) ( - 4^/(i)^ k. (27) 
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Adding Eqs. ([26|) and ([27|) . we see that the terms associated with f{t) cancel and lead to the 
expression 

Lem + Lm = mgRtk. (28) 
as it should, being the change in angular momentum due to the external torque on the system. 

IV. ENERGY BALANCE 

As the hanging mass drops a distance H, the loss in potential energy is 

U = mgH = mgR(t){t), (29) 

where (j){t) is given by Eq. (jl9p . At first sight it may seem that this potential energy is distributed 
into two channels : (i) gain in kinetic energy of system (mass+cylinder) and (ii) electrical energy 
stored in the system. Let us evaluate the contributions of these two channels. 

The total kinetic energy after the hanging mass has dropped a distance H = R(f){t) is 

1 2 1 2 
K = —luj H — mv 

2 2 

= 2 ^ ■ 

where v is the linear speed of the falling mass and u is given by Eq. (|14p . On the other hand the 
stored electrical energy (Ug,) is 

R 



' El 



?7e = ^ / El2^lr'dr' (31) 
^ ./n 



f\t). (32) 



4rM 

Note that the integration is up to R only (See Eq. dS])). However adding Eqs. ([30|) and ([32]) yields 



'2r4r3 



U -{K + Ue) = {I + mR')^^ 2Tj{t)-{f{t)-tY . (33) 



The right hand side is evidently not zero. We need to evaluate the Joule loss. This is given by 



Uj = - I I J''27rlr' dr'dt', (34) 



t rR 
o" Jo JO 

where the current density is given by Eq. (j23p . The heat loss in a conductor is due to the collision 
of the charge carriers with the surrounding atoms. The collisions result in the dissipation of kinetic 
energy of the charge carriers in the form of heat. This motion of charge carriers arises from two 
mutually opposing processes: (i) the electric field Ea given by Eq. ([2T|) and (ii) the Lorentz force. 
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Note that in the evaluation of angular momentum (Eq. (|25p l we considered only electric field (Ea). 
The point to appreciate is that in the case of energy conservation we need to also consider the 
Lorentz force since it results in the motion of charge with consequent loss in kinetic energy due to 
collisions. Substituting the current density from Eq. ([23|) and using Eq. (fT5|) . we obtain 

f,,.^^?!W/V,^(:_,-.7.,V. (35) 

^Tl Jo 

Uj can be evaluated by lengthy and straightforward integration (See Appendix for detailed deriva- 
tion of Uj). It yields precisely the r.h.s. of Eq. (I33p . Thus 

mgH = K + Ue + Uj. 



V. DISCUSSION 

Let us recapitulate the physical processes which are involved. The magnet rotates due to the 
gravitational torque of the falling mass. The free charge in the rotating magnet experiences a 
radially outward force. The resulting radial current is non-solenoidal. This implies that charge 
will accumulate. The electric field due to the charge accumulation will retard the radial motion of 
free charge. The resulting diminished radial current will nevertheless give rise to a counter-torque 
which is of patently electromagnetic origin. Newton's second law can be applied to derive the 
angular velocity and angular acceleration. Next, we demonstrate by detailed bookkeeping that 
energy and angular momentum are properly accounted for. 

One may get a perspective on the problem by obtaining order of magnitude estimates of time 
constants involved. Let us take the radius of the magnet R= 0.1 m, length / = 1 m, density same 
as that of iron (7800 kg-m^^), conductivity a = 10^ S-m^^ and magnetic field i? to be 1 T. Let 
the falling mass be 10 kg. This yields ri ~ 10^^ s, T2 ~ 10""'^'^ s, T3 ~ 10^^ s and Tg ~ 10^"*^^ s. 
Note that t^, = T2. When one uses these numbers then one finds that current density is very 
small and of order 10~^^ A-m~^. Thus the motion of the charges will have a negligible effect on 
the magnetic field. The observable effects are possible if the magnetic field is unusually large as 
happens in magnetars with magnetic field of gigateslas. 10|] However it is important to realise that 



howsoever small the effects due to magnetic field may be they are crucial in accounting for the 
angular momentum and energy balance. 
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Appendix 

We substitute the expression for current density from Eq. ([23]) into Eq. (plj) and get 

Tl Jo Jo 

The integrand associated with the t variable can be simplified with the help of Eq. (jl5p for f{t) 
and this yields 

^^2^2 / 2 4 ,,,,,, 
Uj = ^^^4 ' tI(\ - e-* l^^fdi! 

Expanding and integrating the squared term and using the time constant T3 (Eq. (I13p ) 

f/j = (/ + mi?2)4^ + ?(1 - e-^*/-) - 2r,(l - e"*/^' 
rfrs V 2 

Using Eq. (llSp once again 



= (, + + iziw f 2 - ^ I - 2(t - /(*)) 



TiTs \ 2re 



(/ + mR^-^ \ fit) - ' ) (A.2) 



Equation ()A.2p is precisely the r.h.s. of Eq. ([33p . 
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